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Composite materials, with statistically distributed threshold for breakdown of individual elements, 
are considered. During the failure process of such materials under external stress (load or voltage), 
avalanches consisting of simultaneous rupture of several elements occur, with a distribution D(A) 
of the magnitude A of such avalanches. The distribution is typically a power law -D(A) oc A . 
For the systems we study here, a crossover behavior is seen between two power laws, with a small 
exponent £ in the vicinity of complete breakdown and a larger exponent £ for failures away from 
the breakdown point. We demonstrate this analytically for bundles of many fibers where the load 
is uniformly distributed among the surviving fibers. In this case £ = 3/2 near the breakdown point 
and 5 = 5/2 away from it. The latter is known to be the generic behavior. This crossover is a signal 
of imminent catastrophic failure of the material. Near the breakdown point, avalanche statistics 
show nontrivial finite size scaling. We observe similar crossover behavior in a network of electric 
fuses, and find £ = 2 near the catastrophic failure and £ = 3 away from it. For this fuse model 
power dissipation avalanches show a similar crossover near breakdown. 

PACS numbers: 



I. INTRODUCTION 

Burst avalanches play an important role in character- 
ising the fracture-failure phenomena [H 0- 0- 01 • When 
a weak element in a loaded material fails, the increased 
stress on the remaining elements may cause further fail- 
ures, and thereby give a failure avalanche in which sev- 
eral elements fail simultaneously. With further increase 
in the load new avalanches occur. The statistics of such 
avalanches during the entire failure process explore the 
nature of correlations developed within the system. From 
the experimental point of view, failure avalanches are the 
only measurable quantity during the fracture-failure pro- 
cess of composite materials yiuJQ]. Under quasi-static 
loading the system, with some internal load redistribu- 
tion mechanism, gradually approaches the global failure 
point. Such damage and fracture of materials are of im- 
mense interest due to their economic and human costs. 
Therefore, a fundamental challenge is to find methods 
for providing signals that warn of imminent global fail- 
ure. This is of uttermost importance in, e.g., the diamond 
mining industry where sudden failure of a mine is always 
catastrophic. These mines are under continuous acous- 
tic surveillance, but at present there is no meaningful 
acoustic signature of imminent disaster. The same type 
of question is of course central to earthquake prediction 

0M. 

In this paper we will study crossover behavior of fail- 
ure avalanches in the context of two very different mod- 
els where the system gradually approaches global failure 
through several intermediate failure events. We find that 
if a histogram of the number of elements failing simulta- 



neously is recorded, it follows a power law with an expo- 
nent that crosses over from one value to a very different 
value when the system is close to global failure. This 
crossover is, then, the signature of imminent breakdown. 
The first system studied here is a bundle of many fibers 
HQ with stochastically distributed fiber strengths. This 
model is sufficientl y si mple that an analytic treatment is 
feasible Gil El Hill El El El E3 The second system 
is a fuse model p]], a two-dimensional lattice in which the 
bonds are fuses, i.e., ohmic resistors with stochastically 
distributed threshold values. This model must be ana- 
lyzed numerically. Both models exhibit similar crossovers 
as signal of imminent breakdown. 

The paper is organized as follows. In Section II we 
present numerical evidence for the crossover in the fiber 
bundle model, backed up by analytic derivations. We pay 
particular attention to the burst properties just before 
complete breakdown. Cascading failures in a fuse model 
is the theme of Section III. 



II. THE FIBER BUNDLE MODEL 

A. Numerical evidence 

A bundle of many fibers with stochastically distributed 
fiber strengths, and clamped at both ends, is a much- 
studied model @, H El El El El El El El E3 for 

failure avalanches. In its classical version, a ruptured 
fiber carries no load and the increased stresses caused 
by a failed element are shared equally by all the sur- 
viving fibers. The maximal loads x n that the fibers 
n = 1,2, ...,N are able to carry are picked indepen- 
dently with a probability density p(x) : 



Prob(x < x n < x + dx) — p{x) dx. 



(1) 
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A main result for this model is that under mild restric- 
tions on the fiber strength distribution the expected num- 
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ber D(A) of burst avalanches in which A fibers fail si- 
multaneously is governed by a universal power law ^3] 



D(A) oc A 



-i 



(2) 



for large A, with £ = 5/2. However, we will show that 
when the whole bundle is close to breaking down the 
exponent crosses over to a lower value. Such a complete 
breakdown can be estimated as follows: The force F(x) 
that the bundle is able to withstand when all fibers with 
strengths less than x have ruptured, is proportional to 
the number of surviving fibers times the strength, 



F(x) =Nx Q(x), 



where 



Q(x) = / p(x) dx 



(3) 



(4) 



is the expected fraction of fibers with strengths exceeding 
x. As an example, assume the threshold distribution p(x) 
to be uniform in an interval (xo,x m ), 

n (>A - / ( Xm ~ for x ° - x - Xm 

P[X) ~ otherwise. [b) 



In this case we obtain 



F(x) = N 



•Em %0 



(6) 



which has a maximum at x c = x m /2. In general we 
call F(x) the average force and the value corresponding 
to the maximum of F(x) the critical threshold value x c . 
If F(x) given by (0 were the actual force, the bundle 
would break down when x reaches the value x c . By the 
existence of fluctuations, however, the maximum value 
of the force may actually occur at a slightly different 
(probably higher) value of x. 



results for D(A) for the uniform distribution with xq = 
0.9 x c - For comparison, simulation results with x = 
are shown. In both cases D(A) shows a power law decay, 
apparently with an exponent £ = 3/2 for xq — 0.9 x c in 
contrast to the standard exponent £ = 5/2 for the xq = 
case. 

In Fig. 2 we show that the same two exponents appear 
for a much more concentrated threshold distribution, the 
Weibull distribution. We will in the following explain the 
results as a crossover phenomenon. 
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FIG. 2. The distribution of bursts for the Weibull distribution 
Q(x) = exp(— (x — l) 10 ), where 1 < x < oo. Again results for the 
two cases Xq = 1 (squares) and xq = 1.7 (circles) are displayed 
(x c = 1.72858 for this distribution). The figure is based on 50000 
samples with N = 10 6 fibers and the arrow locates the crossover 
point A c ~ 14.6. 



B. Analytical treatment of the crossover 

For a bundle of many fibers the number of bursts of 
length A is given by [l^l 



D(A) 
D(l) io 2 



<V x = 0.9x c 




FIG. 1. The distribution of bursts for the strength distribution 
l5l with xq = and xq = 0.9x c . The figure is based on 50000 
samples with N = 10 6 fibers. 

We want to study burst avalanches when the weakest 
fiber xq, is close to the critical value x c . In Fig. 1 we show 



g(A) . 
N 

where 



A! 



p(x)r(x)[l — r(x)]^ 1 exp [A r(x)] dx, 

(7) 



r(x) = 1 — 



xp(x) 



1 



Q(x) Q(x) dx 



— \xQ{x)] 



(«) 



From the last expression we see that r(x) vanishes at 
the point x c where the average force expression JHJ is 
maximal. If we have a situation in which the weakest 
fiber has its threshold xq just a little below the critical 
value x c the contribution to the integral in the expres- 
sion JJJ for the burst distribution will come from a small 
neighborhood of x c . Since r(x) vanishes at x c it is small 
here, and we may in this narrow interval approximate the 
A-dependent factors in JjJ as follows 



(l-r) A e Ar = exp[A(ln(l-r) + r)] 
= exp[-A(r 2 /2 + C(r 3 ))] 



! cxp [-Ar(x) 2 /2](9) 



3 



We also have 



r(x) » r'(x c )(x - x c ). 



(10) 



Inserting everything into Eq. 0, we obtain to dominat- 
ing order 



N 



\ A-l -A 



A! 



/ p(x a ) r'(x a )(x - x c )e 



-Ar'( Ic ) 2 (^-^) 2 /2 



dx 



A*- 2 e-*p(x c ) r_ Ar ' (xc) 2 (x _ Xc) 2 /2 

|r'(x c )|AI L 
A A - 2 e- A p(x c ) r _ A/Acl 



A! \r'{x )\ 



(11) 



with 



r'(a; c ) 2 (j: c - x ) 2 ' 



(12) 



By use of the Stirling approximation A! ~ 

A A e _A \/27rA - a reasonable approximation even for 
small A, the burst distribution itTT)) may be written as 



-P(A) 



CA" 5 / 2 (l - e' A / z 



with a nonzero constant 

C=(27r)- 1 / 2 p(x c )/\r'(x c )\ 



(13) 



(14) 



We see from 111 -it that there is a crossover at a burst 
length around A c , so that 



-P(A) 



A -3/2 
A -5/2 



for A < A c 
for A > A c 



(15) 
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FIG. 3. The distribution of bursts for the uniform threshold dis- 
tribution with xo = 0.8Cte c .The figure is based on 50000 sam- 
ples with N = 10 6 fibers. The straight lines represent two different 
power laws, and the arrow locates the crossover point A c ~ 12.5. 

We have thus shown the existence of a crossover from 
the generic asymptotic behavior D oc A~ 5 / 2 to the power 
law D oc A -3 / 2 near criticality, i.e., near global break- 
down. The fact that there may be a different burst dis- 
tribution exponent near breakdown has been noted by 
Sornette (see Ref.l and references therein), and observed 



by Zapperi et al [19] for a fuse model. The crossover 
is a universal phenomenon, independent of the threshold 
distribution p(x) . In addition we have located where the 
crossover takes place. 

For the uniform distribution A c = (1 — xq/x c )~ 2 /2, so 
for xq = 0.9 x c , we have A c = 50. The final asymptotic 
behavior is therefore not visible in Fig. 1. The crossover 
is seen better for xq = 0.8x c , as in Fig. 3. Now a crossover 
is clearly observed near A = A c = 12.5, as expected. 

The simulation results shown in the figures are based 
on averaging over a large number of fiber bundles with 
moderate N. For applications it is important that 
crossover signals are seen also in a single sample. We 
show in Fig. 4 that equally clear power laws are seen for 
a single fiber bundle when N is large. 
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FIG. 4. Avalanche distribution for the uniform threshold distri- 
bution J5j for a single fiber bundle with 10 7 fiber: all avalanches 
(squares) and avalanches near the critical point (circles). Dotted 
straight lines are the best fits to the power laws. 



C. Sampling a finite interval 

Above we have explained the crossover in burst distri- 
butions in which all bursts up to complete breakdown are 
counted. For the purpose of finding signals of imminent 
global failure one must, of course, determine burst distri- 
butions short of complete breakdown. Consequently we 
are interested in sampling a finite interval (xo,Xf), with 
Xf < x c . When the interval is in the neighborhood of x c 
we have, as in Eq.(ll) 



-P(A) 
N 



A A - 2 e- A p{x c ) 



-r'(x c ) 2 (x-x a ) 2 A/2 



x 



(16) 



|r'(z c )|A! 



with a = 2/r'(x c ) 2 . 

This shows a crossover: 



-P(A) 
N 



C A 
CA" 



~ 3 / 2 for A < a/(x c - x ) 2 
5 / 2 for a/(x c - x ) 2 < A < a/(x c - x f ) 2 , 

(18) 
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with a final exponential behavior when A ^> a/ (x c —x f) 2 . 
Here C = Ca' 1 [(x c - x f ) 2 - (x c - x Q ) 2 ] . 

The 3/2 power law will be seen only when the begin- 
ning of the interval, xq, is close enough to the critical 
value x c to create a sizeable range of bursts obeying this 
power law. Observing the 3/2 power law is therefore a 
signal of imminent system breakdown. 



D. Burst avalanches at criticality 

Precisely at criticality (xq = x c ) we have A c = oo, and 
consequently the £ = 5/2 power law is no longer present. 
We will now argue, using a random walk representation, 
that at criticality the burst distribution follows a 3/2 
power law. The load on the bundle when the fcth fiber 
with strength Xk is about to fail is proportional to 



F k =x k (N-k + l). 



(19) 



The expectation value of this is the average force equation 
(jHJ. At criticality the F is, on the average, stationary. It 
is, however, the fluctuations of this load that now deter- 
mines the size of the bursts. It has been shown ^lj that 
the probability p(f) df that the difference Fk+i — Fk lies 
in the interval (/, / + df) is given by 



P(f) 



l-r(xk) -(l-r(x k ))(l+f/xu) for f > 

T. I , •> 



-Xk 



for / < -x k 



(20) 

where r(x) is given by Eq. |H1 At criticality r — 0, result- 
ing in 



Pc(f) 



.-1 p-1 







for / > — x c 
for / < — x c 



(21) 



This can be considered as the step probability in a ran- 
dom walk. The random walk is unsymmetrical, but un- 
biased, (/) =0, as it should be at criticality. 

A first burst of size A corresponds to a random walk 
in which the position after each of the first A — 1 steps 
is lower than the starting point, but after step no. A the 
position of the walker exceeds the starting point. The 
probability of this equals 

Prob(A) = f° x p{h)dh JT/ 1 p{h)dh fZi 1 ~ h P(h)df 3 . . . 

-fl— /2--- — /a- 



p(/a-i)#a-i / K/a)<2/A' 
J — fi— fa—— /a-i 

(22) 



The last integral is easy. By means of 112111 we have 



-fi—h——fA- 
Since p(f)e^ x = e~ x jx we end up with 



(23) 



Prob(A) = e- A dh /J/ 1 dh fir 12 df 



-h-h 



-/l-/2----/A-2 



df 



A-l- 



(24) 



For simplicity we have put x c — 1, since the quantity x c 
simply determines the scale of the steps, and here it is 
only relative step lengths that matters. 

In Appendix A we have evaluated the expression ll22ll , 
with the result 



Prob(A) 



; -a a a-i 
A! 



l - A~ 3 / 2 , (25) 



and also shown that these probabilities satisfy 

oo 

Prob(A) = 1. (26) 
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FIG. 5. Distribution of first bursts (squares) and total bursts 
(circles) for the critical strength distribution J5j with xq = x c . The 
simulation results are based on 10 6 samples with N = 80000 fibers. 
The 'star' symbol stands for the analytic result 1251 . 



The result <2fit is strictly applicable only in the limit 
N — > oo; for finite N the sum has to be slightly less than 
unity. For one thing A cannot exceed N, and breaking off 
the sum at A = N would decrease the sum by an amount 
of order N -1 / 2 . In reality, the sum deviates from unity 
by an amount of order N^ 1 / 3 (see the following section). 

The simulation results in Fig. 5 are in excellent agree- 
ment with the distribution 112511 . At the completion of a 
burst the force, i.e., the excursion of the random walk, 
is larger than all previous values. Therefore one may use 
this point as a new starting point to find, by the same 
calculation, the distribution of the next burst, etc. Con- 
sequently the complete burst distribution is essentially 
oc A" 3 / 2 as expected. In the next section we study the 
burst distribution at criticality in more detail, in partic- 
ular its dependence upon the bundle size N. 



E. 



Finite-size effects at criticality 



When the fiber bundle is sub-critical the average num- 
ber of bursts of a given length will be proportional to 
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the bundle size N. When the bundle is critical this is no 
longer so. Each burst will produce a non-negligible weak- 
ening of the bundle, so that the bundle will be slightly 
more supercritical. Then the probability of a total break- 
down will increase, and the probability of a burst of finite 
length decreases. 

To study this quantitatively we therefore have to spec- 
ify not only the size A of a burst, but also if it is the 
first burst after starting, the second, the third, etc. Let 
P„(A) be the number of bursts of size A that occur as 
the nth burst. If we start precisely at criticality, we have 
already calculated 



^i(A) 



vA-l 



A! 



(27) 



by Ea. 1250 . We will in particular study how the proba- 
bility decreases with n, so we form the ratios 



A(A)- 



(28) 



We start by investigating the A-dependence of these 
ratios, and for simplicity we work with the critical uni- 
form threshold distribution throughout this subsection. 
In Fig. 6 we have plotted R n (A) versus A. The ratios 
H28I1 depend upon n, but surprising enough we cannot 
detect any systematic dependence on A. We may there- 
fore obtain the dependence upon n and N by sticking to 
one fixed A; for simplicity we take A = 1. 
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FIG. 6. Simulation results for the A-dependence of the ratios 
^281 . We have used the uniform strength distribution at criticality. 
The figure is based on 10 7 samples with N = 40000 fibers. The 
straight lines are the average values. 

In Fig. 7 we have plotted i?„(l) for four different values 
of N and for n = 2,3,4 and 5. The figure shows that 



R n (l) for each value of n apparently depends linearly 
on TV -1 / 3 and that 1 — i? n (l) apparently is proportional 
to n — 1. Empirically the data can reasonably well be 
represented by 



Rn(l) = 1 - 1.27 (n- 1)N' 



-1/3 



(29) 



More generally we may assume that the linear function 
is a limiting form of a more general function: 



R n (l) = F(x), 



with x = (n — 1)N 



-1/3 



(30) 



For 



where F(0) — 1 and F(x) ~ 1 — 1.27a; for small x. 
large x we expect F{x) to approach zero. 

If 13011 is correct we should have a data collapse onto 
the single curve F(x). Fig. 8 shows that the data collapse 
works well. In order to test the function F(x) beyond its 
initial linear behavior, we have added a few points with 
larger values of a;. In addition to the results of Fig. 7 
(A) we have obtained results for n = 10, 20 and 30, with 
N = 5000, N = 10000 and N = 80000 (Fig. 7 (B)). 
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FIG. 7. Simulation results for the ratios -Rn(l) as function of 
TV -1 / 3 (A) and as function of n (B). The dotted lines represent the 
functional forms 1291 (in A) and ^311 (in B). The results are based 
on 10 7 samples, except for N = 80000 (10 6 samples). 
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F{x) is seen to decrease towards zero for increasing x, 
and the empirical expression 



F(x) 



„2.54x 



(31) 



seems to represent the data in Fig. 8 very well. 

The relation I-IOII implies some interesting conse- 
quences. Let us first consider the total number of burst 
of a given size: 

D(A) = ^P„(A) = P 1 (A)^fl„(A)~P 1 (A)^ J R„(l) 

= Pi(A)J2 F l( n - l)N~ 1/3 ] ~ Pi (A) j °° F[(n - 1)N~ 1/3 ] dn 

= P 1 (A)N 1/3 f°° F(x) dx. (32) 
J o 



We have used that i? n (A) is essentially independent of 
A (Fig. 6), and that due to smallness of N' 1 ^ 3 we may 
replace summation by integration. 
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FIG. 8. The data collapse onto a single curve F{x) (dot- 
ted line represents Eq. 28) where x = (n — 1) x iV -1 / 3 . We 
have taken N = 5000, 10000, 20000, 40000, 80000 and n = 
2, 3, 4, 5, 10, 20, 30. Averages are taken over 10 7 samples, ex- 
cept for N = 80000 (10 6 samples). 

The conclusion is that the total number of bursts 
should scale as TV 1 / 3 . The simulation results presented 
in Table 1 are in excellent agreement with this N- 
dependence. 



Table 1 



N 


5000 


10000 


20000 


40000 


80000 


N-^ 3 D(1)/P(1) 


0.545 


0.546 


0.546 


0.550 


0.550 


N-y a D(2)/P(2) 


0.546 


0.549 


0.550 


0.547 


0.552 


N- 1 ^D(3)/P(3) 


0.542 


0.550 


0.548 


0.547 


0.546 


iV- 1 / 3 D(4)/P(4) 


0.538 


0.548 


0.544 


0.548 


0.546 


iV- 1 / 3 L>(5)/P(5) 


0.542 


0.545 


0.546 


0.548 


0.550 


jV^ 3 D(failure) 


1.2154 


1.2319 


1.2346 


1.2293 


1.2307 



According to (1321) the numbers in the table should all 
be the same (or nearly the same) and represent the in- 
tegral of F. The integral of the empirical representation 



113 111 of F(x) equals 0.546, in close agreement with the 
results in the table. 

We have also recorded the number D (failure) of im- 
mediate failures of the fiber bundle (i.e., with no finite 
bursts at all). We have presented the numbers in Ta- 
ble 1. The number of immediate failures decreases with 
increasing N as iV" 1 / 3 . The reason for the decrease is 
that in a large bundle it is more probable to find a fiber 
sufficiently strong to prevent immediate failure. 

It remains a challenge to derive these finite-size scaling 
results analytically. 



III. 



BURST AVALANCHES IN THE FUSE 
MODEL 




FIG. 9. Burst distribution in the fuse model: System size is 
100 x 100 and averages are taken for 300 samples . On the aver- 
age, catastrophic failure sets in after 2097 fuses have blown. The 
circles denote the burst distribution measured throughout the en- 
tire breakdown process. The squares denote the burst distribution 
based on bursts appearing after the first 1000 fuses have blown. 
The triangles denote the burst distribution after 2090 fuses have 
blown. The two straight lines indicate power laws with exponents 
£ = 3 and 5 = 2, respectively. 

Let us test the crossover phenomenon in a more com- 
plex situation than for fiber bundles. We have studied 
burst distributions in the fuse model It consists of 
a lattice in which each bond is a fuse, i.e., an ohmic re- 
sistor as long as the electric current it carries is below 
a threshold value. If the threshold is passed, the fuse 
burns out irreversibly. The threshold t of each bond is 
drawn from an uncorrelated distribution p(t) . The lattice 
is placed between electrical bus bars and an increasing 
current is passed through it. Numerically, the Kirchhoff 
equations are solved with a voltage difference between 
the bus bars set to unity. The ratio between current ij 
and threshold tj for each bond j is calculated and the 
bond having the largest value, maxj(ij /tj) is identified 
and subsequently irreversibly removed. The lattice is a 
two-dimensional square one placed at 45° with regards 
to the bus bars. The threshold distribution is uniform 
on the unit interval. All fuses have the same resistance. 



7 



The burst distribution follows the power law 101 with 
£ = 3, which is consistent with the value reported in 
recent studies 0,0. We show the histogram in Fig. 
9. With a system size of 100 x 100, 2097 fuses blow 
on the average before catastrophic failure sets in. When 
measuring the burst distribution only after the first 2090 
fuses have blown, a different power law is found, this time 
with £ = 2. After 1000 blown fuses, on the other hand, 
£ remains the same as for the histogram recording the 
entire failure process, see Fig. 9. Zapperi et al [19], who 
study the fuse model on the diamond and the triangular 
lattices, find significant variation with the lattice type. 
Their exponent values for the diamond lattice are 2.75 
and 1.90, not very different from the values 3.0 and 2.0 
in our Fig. 9. 

In Fig. 10, we show the power dissipation E in the 
network as a function of the number of blown fuses. The 
dissipation is given as the product of the voltage drop 
across the network V times the total current that flows 
through it. In Fig. 11, we show the power dissipation as 
a function of the total current. The breakdown process 
starts by following the lower curve, and follows the upper 
curve returning to the origin. It is interesting to note the 
linearity of the unstable branch of this curve. In Fig. 12, 
we record the avalanche distribution for power dissipa- 
tion, Dd(A). Recording, as before, the avalanche distri- 
bution throughout the entire process and recording only 
close to the point at which the system catastrophically 
fails, result in two power laws, with exponents £ = 2.7 
and £ = 1.9, respectively. It is interesting to note that 
in this case there is not a difference of unity between the 
two exponents. The power dissipation in the fuse model 
corresponds to the stored elastic energy in a network of 
elastic elements. Hence, the power dissipation avalanche 
histogram would in the mechanical system correspond to 
the released energy. Such a mechanical system would 
serve as a simple model for earthquakes. 
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The Gutenberg-Richter law 0, 0- Q| relating the fre- 
quency of earthquakes with their magnitude is essentially 
a measure of the elastic energy released in the earth's 
crust, as the magnitude of an earthquake is the loga- 
rithm of the elastic energy released. Hence, the power 
dissipation avalanche histogram Aj(A) in the fuse model 
corresponds to the quantity that the Gutenberg-Richter 
law addresses in seismology. Furthermore, the power law 
character of Dd(A) is consistent with the form of the 
Gutenberg-Richter law. It is then intriguing that there 
is a change in exponent £ also for this quantity when 
failure is imminent. 
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FIG. 11. Power dissipation E as a function of the total 
current I flowing in the fuse model. The system size and 
number of samples are the same as in Fig. 9. 
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FIG. 10. Power dissipation E as a function of the number of 
broken bonds in the fuse model. The system size and number 
of samples are the same as in Fig. 9. 



FIG. 12. The power dissipation avalanche histogram Dd{A) 
in the fuse model. The slopes of the two straight lines are 
—2.7 and —1.9, respectively. The circles show the histogram 
of avalanches recorded through the entire process, whereas 
the squares show the histogram recorded only after 2090 fuses 
have blown. The system size and number of samples are the 
same as in Fig. 9. 



8 



IV. CONCLUDING REMARKS 



we have 



We have studied the avalanche distribution D(A) oc 
A"^ in the fiber bundle model, and have shown ana- 
lytically that close to complete breakdown it exhibits a 
crossover behavior between two power laws with expo- 
nents £ = 5/2 and £ = 3/2. This crossover behavior 
is universal in the sense that, under mild assumptions, 
it does not depend on the statistical distribution of the 
thresholds. In the critical situation an argument based on 
a unbiased unsymmetrical random- walk scenario explains 
the exponent £ = 3/2. Near criticality the avalanche dis- 
tribution depends on the system size in a nontrivial way. 
For this case we present quantitative results that may be 
summarized by a finite-size scaling function [Eq. l3lj . 

The crossover behavior is not limited to the fiber bun- 
dle model. We show numerically that the same crossover 
phenomenon occurs in the two-dimensional fuse model. 
The exponents are different, though, £ = 2 near break- 
down and £ = 3 away from it. For this fuse model 
the power dissipation avalanches show a crossover, with 
power law exponents £ = 2.7 and £ = 1.9. Such crossovers 
signal that catastrophic failure is imminent, and has 
therefore a strong potential as a useful detection tool. 
Some of the present results have already been published w ith 
as a letter (l7l | 



Prob(A) =e- A V A -i(0). 



(38) 



Equation l|37|) can be solved by iteration. By calculating 
the first polynomials V n (y) one is led to assume 



d-i 



(i + l)y\ (39) 



Suppose this is valid up to some value of d— 1. Then use 
(53 to compute Vd- The integration is trivial, leaving 



1 d-l 

w = * E 



d- 1 



jd—i—l 



(1+2/)' 



d 

' ^ \ i 1 

i=0 v 7 

, d-l i+1 , , 

l)- W'' .')('" ) (/ ./ , i 

^" i=0 m=0 ^ 

d d—1 / , / , 

1 ^-v \ ^ / d—1 \ / i + 1 



d-l \ (i + 1 
i ) \ m 



d\ 



m— i— m — 1 

\ E S(m) 



jd—i — l , 



m=0 
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APPENDIX A: PROOF OF Ea.(f25|l 

We evaluate here the multiple integral in Eq. 1(221) 

b(A) = e" A J d/i J df 2 J i dh...J 

We introduce the new variables 



S(m)= ]T 



d-l \ ( i + 1 \ ,,/- , i 



Since 



i + 1 
rn 



m — 1 



and 



Pro 



(33) 



2/1 = -/l 

2/2 = -fl - /2 

2/A-l = ~fl - /2 



we may write 

d-l 

S(m) = 

i=0 



= 



d-l 



for b < or 6 > 



(40) 



(41) 



(42) 



(43) 



m — 1 



jd— i— 1 



(44) 

/a-i, (34) To evaluate S'(to) we differentiate the binomial expres- 



satisfying 



< yi < 1 

< 2/2 < 1 + 2/i 

0< Ui <l + y i+ i £ = 2,3, . . . , A - 1. (35) 

Then 

Prob(A) = e a d Vl dy 2 ... dy A -l- (36) 

./ ■/ ■/ 

Defining Vo(2/) = 1, and 

ri+y 

V d (y)= / V d -i(z)dz, (37) 
Jo 



E( d 7 1 ) :cl = (1+ ^ 1 



(45) 



to times with respect to x: 

d-l , . 

J2 ( 1 ) ■ • ■ (t-m+1) a;* - "* = (d-l)(d-2) . . . (d-m)(l+x) d_1 " 

(46) 

or 

rf-i 

E 

i=0 



d-l 



(d-l) 



(d-m-1)! 



(1+x) 



d-l- 1 



(47) 
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Putting now x — 1/d, and multiplying both sides of the which is Ea . (12511 in the main text. 



equation by d d 1 m /m\, we obtain 

d-l 



s ( m ) = 



d-l 



d 



(d-l)! 



(d — m — l)!m! 



(l+d)" 



By Eq.gU then 

S(m) = s(m)+s(m—l) 
Finally Eq.JUjJ gives 



(48) 



(l+d) d - m (l+m)/(l+d). 

(49) 



(l + d)! 



E 



(d+i) d -" l (i+Hy m - (so) 



Since this is in accordance with the assumption (|39)) , and 
since iKfflU is correct for d = 2, the induction proof works. 
Finally, using l(38jl . the probability we seek is 



Prob(A) = e- A V A -i(0) = 



-A A A-i 



Let us also sum these probabilities over all burst 
lengths, 



A! 



T — 



dS A - 



(52) 



We may now appeal directly to the theorem of Lagrange 
20] to conclude that the sum satisfies the equation 



S = e-V. 



(53) 



Since S e s is always less or equal to e 1 for nonnegative 
S. we must have 



5 = 1, 



(54) 



(51) which is En. ( 12611 in the main text. 
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